Bulk superconducting phase with a full energy gap in the doped topological insulator 

Cu x .Bi 2 Se 3 
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The superconductivity recently found in the doped topological insulator CuzEfeSes offers a great 
opportunity to search for a topological superconductor. We have successfully prepared a single- 
crystal sample with a large shielding fraction and measured the specific-heat anomaly associated 
with the superconductivity. The temperature dependence of the specific heat suggests a fully- 
gapped, strong-coupling superconducting state, but the BCS theory is not in full agreement with 
the data, which hints at a possible unconventional pairing in CuzBi2Se3. Also, the evaluated effective 
mass of 2.6m e (m e is the free electron mass) points to a large mass enhancement in this material. 

PACS numbers: 74.25. Bt; 74.70.Ad; 74.62.Dh; 74.25. Op 



In the past two years, the three-dimensional (3D) topo- 
logical insulator (TI) is attracting a lot of interest as a 
new state of matter ■ It is characterized by the ex- 
istence of a gaplcss surface state that emerges because 
of the non-trivial Z2 topology of the insulating bulk 
state and is protected against backscattering by time- 
reversal symmetry The discovery of the 3D TI stimu- 
lated the search for a superconducting (SC) analogue, 
a time-reversal-invariant topological superconductor 0- 
[§], which is characterized by a fully-gapped, odd- parity 
pairing state that leads to the emergence of gapless Ma- 
jorana surface states. Such a SC phase has implications 
on topological quantum computing [i, T3-12| because of 
the non-Abelian Majorana bound state expected to ap- 
pear in the vortex core. However, a concrete example of 
such a topological superconductor is currently unknown. 

In this context, the superconductivity recently found 
[lit in Cuj;Bi2Se3 is very interesting. Bi2Se3 is a "second- 
generation" TI that has a relatively large (~0. 3 eV) band 



gap and a simple surface-state structure 13|,LL|. Surpris- 
ingly, when Cu is introduced to this system with the nom- 
inal formula Cu^I^Ses, superconductivity with a maxi- 
mum transition temperature T c of 3.8 K was observed for 
the doping range 0.10 < x < 0.30, even though the bulk 
carrier density n was only ~ 10 20 cm -3 Note that 
this T c is uncharacteristically high for such a low n [lj|. 
Furthermore, the topological surface state was found to 
be well-separated from the doped bulk conduction band 
in Cua;Bi2Se3 [T3] - So far, apparent SC shielding fractions 
of only up to 20% have been achieved and the resistivity 
always remained finite [IB , 17 1 , leaving some doubt about 
the bulk nature of the superconductivity in this system. 
Nonetheless, if this superconductivity is indeed a bulk 
property of carrier-doped Bi 2 Se3, it has a profound im- 
plication on the search of topological superconductors, 
because (i) it is a potential candidate [HI to realize a 3D 
topological superconductor, and (ii) if its bulk turns out 
to be an ordinary s-wave superconductor, the topological 
surface state may turn into a 2D topological supercon- 
ductor as a result of a SC proximity effect Therefore, 



it is important to confirm whether the superconductivity 
is really occurring in the bulk of Cu x Bi 2 Se3 and, if so, to 
elucidate the fundamental nature of its SC state. 

Bi 2 Sc 3 has a layered crystal structure (R3m, space 
group 166) consisting of stacked Se-Bi-Se-Bi-Se quintu- 
ples that are only weakly van-der-Waals bonded to each 
other. We call the rhombohcdral [111] direction the c 
axis and the (111) plane the ah plane. When Cu is intro- 
duced into Bi 2 Sc3, it may either intercalate as Cu 1+ into 
the van-der-Waals gaps and act as a donor, or replace Bi 
as a substitutional impurity and act as an acceptor [l9j |: 
hence, Cu is an ambipolar dopant 2^, 21 1. The nominal 
formula of Cttj;Bi2Se3 suggests that most Cu atoms in 
this SC material occupy the intercalation sites; however, 



the reported carrier density of ~10 cm" J [15( corre- 
sponds to only ~1% of electron doping, which is much 
smaller than that expected from the x value. This dis- 
crepancy suggests either that most of the intercalated 
Cu ions remain inactive as donors, or that substitution 
of Bi with Cu also occurs in this material and it almost 
compensates the electrons doped by the intercalated Cu. 
Partly related to such an uncontrollability of the Cu 
atoms in Cua;Bi2Se3, the quality of the SC samples has 
been poor as mentioned above, and improvements in the 
sample quality are indispensable for a solid understand- 
ing of the SC state in this material. 

In this Letter, we report a comprehensive study of the 
basic SC properties of a Cu-intercalated Bi2Se3 single 
crystal by means of resistivity, magnetization, and spe- 
cific heat measurements. For the first time in this mate- 
rial, we observed zero-resistivity and a specific-heat jump 
at the SC transition. The apparent shielding fraction of 
our sample exceeds 40%, and the specific- heat data con- 
firms the bulk nature of the superconductivity. Most 
importantly, the temperature dependence of the specific 
heat suggests a fully-gapped, strong-coupling SC state, 
but the data do not fully agree with the strong-coupling 
BCS calculation. This suggests that the pairing symme- 
try may not be simple isotropic s-wavc. Furthermore, the 
effective mass is found to be 2.6to c (m e is the free electron 



2 



mass), suggesting a change in the bulk band curvature. 

Single crystals of Bi2Se3 were grown by melting stoi- 
chiometric amounts of elemental shots of Bi (99.9999%) 
and Se (99.999%) in sealed evacuated quartz glass tubes 
at 800°C for 48 h, followed by a slow cooling to 550°C 
over 48 h and keeping at that temperature for 24 h. The 
crystals were cleaved and cut into rectangular pieces, and 
then the Cu intercalation was done by an electrochemi- 
cal technique under inert atmosphere inside a glove box, 
using Cul reagent in CH3CN solvent. The sample was 
wound with a Cu wire, and a Cu stick was used as counter 
and reference electrode. The current was fixed at typi- 
cally 10 fiA. The concentration of intercalated Cu was de- 
termined from the weight change before and after the in- 
tercalation process, and the sample was briefly annealed 
afterward. The particular sample used in the present 
study was 3.9x1.6 mm 2 in the ab plane with a thick- 
ness of 0.40 mm, and its Cu concentration was x = 0.29. 
In fact, by employing electrochemical intercalation, we 
found that samples with x of up to ~0.5 become super- 
conducting; the precise phase diagram is currently under 
investigation. 

The magnetization data were measured with a 
commercial SQUID magnetometer (Quantum Design, 
MPMS). The resistivity p xx and the Hall resistivity p yx 
were measured by a standard six-probe technique, where 
the electrical current was applied in the ab plane. The 
specific heat c p data were taken by a relaxation-time 
method using a commercial system (Quantum Design, 
PPMS); the addenda signal was measured before mount- 
ing the sample and was duly subtracted from the mea- 
sured signal. The c p measurement was done in zero-field 
(for the SC state) and in 2 T applied along the c axis (for 
the normal state), and the change of the addenda signal 
between the two were found to be negligible. 

Figure 1 shows the results of the transport measure- 
ments. In zero field, the onset of the SC transition oc- 
curs at 3.6 K and zero-resistivity was observed at 2.8 K 
[Fig. 1(a)]. From the magnetic- field (B) dependence of 
p xx [Figs. 1(c) and (d)], the upper critical field B c .y(T) 
defined as the midpoint of the resistivity transition [22| 
is obtained for the two principal field directions [Fig. 
1 (e)] , and the extrapolation to K using the Werthamer- 
Helfand-Hohenbcrg (WHH) theory gives i? C 2,±(0) and 
-Bc2,||(0) (perpendicular and parallel to the ab plane) of 
1.71 and 3.02 T, respectively [Fig. 1(c)]. As shown in the 
inset of Fig. 1(b), p yx is completely linear in B, suggest- 
ing the dominance of only one type of bulk carriers. The 
Hall coefficient Rh was found to be only weakly tempera- 
ture dependent [Fig. 1(b)] and gives the electron density 
n = 1.3xl0 20 cm~ 3 . 

Figure 2 summarizes the magnetization M measure- 
ments. To minimize the effect of the demagnetization 
factor, those measurements were made for B \\ ab [see 
the sketch in Fig. 2(a)]. The temperature dependence 
of the diamagnetic shielding fraction for the zero-ficld- 
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FIG. 1. (color online) (a) p xx (T) data of the Cuo.29Bi2Se3 
sample, (b) Temperature dependence of Rh; inset shows the 
Pyx(B) data at 5 K. (c,d) p xx (B) data for B \\ ab and B _L ab, 
respectively, (e) B C 2 vs. T phase diagram determined from 
the midpoint in p xx (B) at various temperatures; dashed lines 
show the WHH behavior. The midpoint definition for B C 2 
gives T c = 3.2 K consistent with the M(T) data. 
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FIG. 2. (color online) (a) Temperature dependence of the ap- 
parent shielding fraction of Cuo.29Bi2Se3 measured in B — 0.2 
mT || ab. (b) M(B) curves at 1.8 and 2.4 K after subtracting 
the diamagnetic background, (c) Initial M(B) behavior af- 
ter ZFC to various temperatures. Arrows mark the position 
of S cl j|. Note the very small magnetic-field scale, (d) Plots 
of AM = M — aB, where a is the initial slope, and the de- 
termination of -B c i. || shown by arrows, (e) B cl 11 vs. T phase 
diagram; the solid line is a fit within the local dirty limit. 
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FIG. 3. (color online) (a) c p (T)/T data measured in and 2 
T applied along the c axis; the latter represent the normal- 
state behavior. The dashed line is a fit to the 2-T data using 
the standard Debye formula, (b) Electronic term c e i/T in 
T obtained after subtracting the phonon term determined in 
2 T. The dotted line shows the calculated c a \/T curve given 
by strong-coupling BCS theory with a — 1.9; the dashed line 
is the BCS curve for a = 2.3, which is obtained from B c , No, 
and T c . The horizontal dash-dotted line denotes the value of 
7 n , and its breakdown to 7 S and 7 rcs is indicated. 



cooled (ZFC) and field-cooled (FC) measurements are 
shown in Fig. 2(a), where the onset of the Meissner signal 
occurs at T c = 3.2 K and the apparent shielding fraction 
reaches 43% at 1.8 K 



23|. Note that this Meissner T c 



corresponds to the midpoint of the resistivity transition. 
Neither ZFC nor FC data saturate at 1.8 K. 

Magnetization M(B) curves are shown in Figs. 2(b) 
and (c); each data set was obtained after cooling to its 
respective temperature from above T c in zero field, and 
the background diamagnetism, which can be easily de- 
termined at B > -B c2 .||, is subtracted from the data. As 
already noted by Hor et al. [lEj . the lower critical field 
B c i is very small: Using the deviation of the M(B) curve 
from its initial linear behavior as a measure of -B c i,|| [Fig. 
2(d)], we obtained the B c i^(T) data shown in Fig. 2(e). 
To determine the 0-K limit, we used B c i oc 1/A 2 ff oc 
f[A(T)/A(0)) tanh(A(T)/2fc B T)] for the local dirty limit 
[25j to fit the extracted data points (A c ff is the effective 
penetration depth and A is the SC gap [i^]), and ob- 
tained B*^%, (0) = 0.43 mT. For the quantitative analysis 
discussed later, this apparent value was corrected for the 
demagnetization effect, though it is small for B \\ ab: Us- 



ing the approximation given for the slab geometry [24|, 
we obtain S cl .|| (0) = B^ p « (0) /tanh^/0.366/a = 0.45 mT, 
where b/a = 3.9/0.40 in our case. Note that the flux pin- 
ning in the present system is weak as evidenced by the 
low irreversibility field of -0.1 T at 1.8 K [Fig. 2(b)]. 

The temperature dependence of c p is shown in Fig. 
3(a) as c p /T vs. T for the SC state (B = T) and 
the normal state achieved by applying B _L ab of 2 T 
(> B C 2.±). As shown by the dotted line in Fig. 3(a), 
a conventional Debye fit to the normal-state data below 



4 K using c p = c el + c ph = 7n T + A 3 T 3 + A 5 T 5 , with 
the normal-state specific-heat coefficient 7n and the co- 
efficients of the phononic contribution A% and A$, yields 
a good description of the data. The obtained parame- 
ters are 7„ =1.95 mJ/molK 2 , A 3 = 2.22 mJ/molK 4 [13], 
and A$ = 0.05 mJ/molK 6 . Subtracting the phononic 
contribution from the zero-field data gives the electronic 
specific heat c e i in the SC state plotted in Fig. 3(b), re- 
vealing a clear jump around T c . This provides compelling 
evidence for bulk superconductivity in CuzI^Sea. In 
passing, we note that our c p data in 2 T do not exhibit 
any Schottky anomaly related to electron spins, suggest- 
ing that there is no local moment possibly associated with 
Cu 2+ ions. 

From the above results, one can estimate various basic 
parameters. Assuming a single spherical Fermi surface, 
one obtains the Fermi wave number fcp = (S^n) 1 ^ 3 = 
1.6 nm -1 . The effective mass m* is evaluated as m* = 
(3fi. 2 7n )/(U mo ifc 2 5 A:F) = 2.6m c , with the molar volume of 
Bi2Se3 Vmoi w 85 cm 3 /mol. Note that the effective mass 
of pristine Bi2Se3 is — 0.2m e [2^], so there is an order-of- 
magnitude mass enhancement in Cu^E^Sea [29j |. Since 
electron correlations are weak in Bi2Se3, the origin of 
this enhancement is most likely a change in the band 
curvature near the Fermi level. From -B C 2,± = 1-71 T, 
the coherence length £ a fc = \/^o/{2ttB c2 ,±) = 13.9 nm is 
obtained, while from B c2 ^ we use £ a &£c = $0 /(27T-B C 2, \\ ) 
and obtain £ c = 7.9 nm. Since we have the B c \ value 
only for B \\ ab, we define the effective GL parameter 
K a b = ^Xab^c/Zabtic and use B cl i| = $ ln/t ab /(47rA ab A c ) 
together with B c2 .\\/ B cl ^ ~ 2/i 2 ^/ In K a f, [30} to obtain 
n a f, « 128. We then obtain the thermodynamic critical 
field B c = y/B cl ^B c2 ,\\/ln.K ab = 16.7 mT. 

To analyze c c \/T in the SC state shown in Fig. 3(b), 
we tried to fit the BCS-type temperature dependence to 
the data. Since the simple weak-coupling BCS model 
does not describe the c e \/T data (not shown), we use 
the modified BCS model applicable to strong-coupling 
superconductors as proposed in Ref. [3l[, where it is 
called "a model" with a = Aq/T c and Ao is the SC gap 
size at K. We note that strong coupling means a > 
a BCS = 1-764, and that this model still assumes a fully- 
gapped isotropic s-wave pairing. Using the theoretical 
curve c^j CS of the a model [3l[ , we tried to reproduce the 



experimental data with c c \(T)/T 



It, 



(T)/T. 



Note that the parameter 7rcs is necessary for describing 
the contribution of the non-SC part of the sample j32| : 
also, the theoretical term c^ cs /T is set to yield 7 S (= 
7n - 7rcs) at T > T c . 

It turned out that with a = 1.9, 7rC s = 0.6 mJ/molK 2 , 
and 7 S = 1.185 mJ/molK 2 , the experimental data is rea- 
sonably well reproduced and the entropy balance is sat- 
isfied, as shown in Fig. 3(b) by the dotted line and the 
dash-dotted line [33|. This result strongly suggests that 
the SC state of Cu K Bi2Sc3 is fully gapped. The resulting 
7n (= 7rcs + 7s) value of 1.785 mJ/molK 2 slightly devi- 
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ates from the 7 n value estimated from the Dcbyc fit to 
the normal-state data in 2 T, 1.95 mJ/molK 2 . This slight 
difference (~9%) might be the result of a possible field 
dependence of the normal-state Sommerfeld parameter, 
which has to be clarified in future studies. 

To gain further insight into the nature of the SC state 
in Cua;Bi2Se3, we examine the implication of the obtained 
7 S value: The density of states (DOS) N is calculated 
from 7 S via N = 7 s /(7r 2 fc|/3) = 1.51 states/eV per unit 
cell, which is large for a low-carrier-density system and 
is in accord with the "high" T c . This value allows us to 
calculate Ao through the expression for the SC conden- 
sation energy ^N A^ = (l/2/i )B 2 . With B c re 16.7 mT 
already calculated, we obtain Ao = 7.3 K which gives 
the coupling strength a = A /T c = 2.3. This exceeds 
the BCS value of 1.764 and hence Cu x Bi2Se3 is a strong- 
coupling superconductor, as was already inferred in our 
analysis of the c e \/T data. More importantly, the a value 
of 2.3 obtained from 7 S is too large to explain the c e i(T) 
data within the strong-coupling BCS theory: As shown in 
Fig. 3(b) with the dashed line, the expected BCS curve 
for a = 2.3 does not agree with the data at all. This 
probably means that the actual temperature dependence 
of A in Cua;Bi2Se3 is different from that of the BCS the- 
ory, which suggests that the pairing symmetry may not 
be the simple isotropic s-wave. Obviously, a direct mea- 
surement of Ao and A(T) is strongly called for. On the 
other hand, the low-temperature behavior of c e \(T) ro- 
bustly indicates the absence of nodes and points to a 
fully-gapped state. It will be interesting to see if the fully- 
gapped, time- reversal- invariant p- wave state proposed for 
Cu a; Bi2Se3 (l8j would provide a satisfactory explanation 
of our data. 

In summary, we report a comprehensive study of the 
superconductivity in Cu T Bi2Se3 by means of resistivity, 
magnetization, and specific-heat measurements on a sin- 
gle crystal with x = 0.29 that shows, for the first time 
in this material, zero-resistivity and a shielding fraction 
of more than 40%. An analysis in the framework of a 
generalized BCS theory leads to the conclusion that the 
superconductivity in this system is fully gapped with a 
possibly non-BCS character. The fully-gapped nature 
qualifies this system as a candidate for a topological su- 
perconductor: Since this system hosts a topological sur- 
face state above T c [17 1 , depending on whether the parity 
of the bulk SC state is even or odd, either the surface or 
the bulk should realize the topological SC state associ- 
ated with intriguing Majorana edge states. 
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